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Abstract 

Unitary space-time modulation using multiple antennas promises reliable communi- 
cation at high transmission rates. The basic principles are well understood and certain 
criteria for designing good unitary constellations have been presented. 

There exist two important design criteria for unitary space time codes. In the 
situation where the signal to noise ratio is large it is well known that the diversity 
product (DP) of a constellation should be as large as possible. It is less known that 
the diversity sum (DS) is a very important design criterion for codes working in a 
low SNR environment. For some special situations, it will be more practical and 
reasonable to consider a constellation optimized at a certain SNR interval. For this 
reason we introduce the diversity function as a general design criterion. So far, no 
general method to design good-performing constellations with large diversity for any 
number of transmit antennas and any transmission rate exists. 

In this paper we propose constellations with suitable structure which allow one 
to construct codes with excellent diversity using geometrical symmetry and numeri- 
cal methods. We also demonstrate how these structured constellations out-perform 
currently existing constellations and explain why the proposed constellation structure 
admit simple decoding algorithm: sphere decoding. The presented design methods 
work for any dimensional constellation and for any transmission rate. Moreover codes 
based on the proposed structure are very flexible and can be optimized for any signal 
to noise ratio. 
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A preliminary version of this paper was presented at 40-th Allerton Conference on Communication, Control, 
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1 Introduction and Model 



One way to acquire reliable transmission with high transmission rate on a wireless channel is 
to use multiple transmit or receive antennas. Either because of rapid changes in the channel 
parameters or because of limited system resources, it is reasonable to assume that both the 
transmitter and the receiver don't know about the channel state information (CSI), i.e. the 
channel is non-coherent. 

In jTH] , Hochwald and Marzetta study unitary space-time modulation. Consider a wireless 
communication system with M transmit antennas and N receive antennas operating in a 
Rayleigh flat-fading channel. We assume time is discrete and at each time slot, signals are 
transmitted simultaneously from the M transmitter antennas. We can further assume that 
the wireless channel is quasi-static over a time block of length T. 

A signal constellation V := {$1, ...,$£,} consists of L matrices having size T x M and 
satisfying T > M and &* k &k = Hi- The last equation simply states that the columns of 
$fc form a "unitary frame", i.e. the column vectors all have unit length in the complex 
vector space C T and the vectors are pairwise orthogonal. The scaled matrices \/T<&k, k = 
1,2, ••• , L, represent the code words used during the transmission. It is known that the 
transmission rate is determined by L and T: 

log 2 (L) 
T 

Let p represent the expected signal-to-noise ratio (SNR) at each receive antenna. The 
basic equation between the received signal R and the transmitted signal \/T<& is given 
through: 

R = 

where the M x N matrix H accounts for the multiplicative complex Gaussian fading coeffi- 
cients and the T x N matrix W accounts for the additive white Gaussian noise. The entries 
h m ,n of the matrix H as well as the entries w tyn of the matrix W are assumed to have a 
statistically independent normal distribution CAf(0, 1). In particular it is assumed that the 
receiver does not know the exact values of either the entries of H or W (other than their 
statistical distribution). 

The decoding task asks for the computation of the most likely sent code word $ given 
the received signal R. Denote by | \\p the Frobenius norm of a matrix. If A = (ajj) then 




the Frobenius norm is defined through = yYlij \ a i,j\ 2 - Under the assumption of the 

above model the maximum likelihood (ML) decoder will have to compute: 

§ml = arg max 

for each received signal R. (See [T9]). 

Let 5 m be the m-th singular value of It has been shown in ^Hj that the 

pairwise probability of mistaking $z for <Jty using maximum likelihood decoding satisfies: 
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-P*;,*,/ = Prob ( choose 3ty | transmitted ) (p) 
= Prob ( choose $2 | $£/ transmitted ) (p) 
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It is a basic design objective to construct constellations V = {<&!,...,<&£,} such that the 
pairwise probabilities small as possible. Mathematically we are dealing with 

an optimization problem with unitary constraints: 

Minimize maxP$ h $ i( with the constraints = I where % = 1, 2, • • • , L. 

Formula (jl.2j) is sometimes referred to as "Chernoff's bound". This formula is easy to 
work with, the exact formula (|l.lj) is in general not easy to work with, although it could be 
useful in the numerical search of good constellations as well. Researchers have been searching 
for constructions where the maximal pairwise probability of -P$ ; ,$,, is as small as possible. 
Of course the pairwise probabilities depend on the chosen signal to noise ratio p and the 
construction of constellations has therefore to be optimized for particular values of the SNR. 

The design objective is slightly simplified if one assumes that transmission operates at 
high SNR situations. In [T£], a design criterion for high SNR is presented and the problem 
has been converted to the design of a finite set of unitary matrices whose diversity product 
is as large as possible. In this special situation several researchers 123 EE EZj came up 
with algebraic constructions and we will say more about this in the next section. 

The main purpose of this paper is to present structured constellation and to develop 
geometrical and numerical procedures which allow one to construct unitary constellations 
with excellent diversity for any set of parameters M, N, T, L and for any signal to noise ratio 
p. The paper is structured as follows. In Section El we introduce the diversity function of a 
constellation. This function depends on the signal to noise ratio and it gives for each value 
p an indication how well the constellation V will perform. For large values of p the diversity 
function is governed by the diversity product, for small values of p it is governed by the 
diversity sum. These concepts are introduced in Section El as well. The introduced concepts 
are illustrated on some well known constellations previously studied in the literature. 

In Section El we first show that randomly constructed codes are fully diverse with proba- 
bility one. Then we start the main task of this paper, namely to parameterize constellations 
which will be efficient for numerical search algorithms. For this purpose we introduce the 
concept of a weak group structure and we classify all weak group structures whose elements 
are normal and positive. 

In Section 13] we investigate an algebraic structure which led to some of the best constel- 
lations which we were able to derive. We also show that in the good-performing codes the 
distance spectrum profile for both the diversity sum and the diversity product are important. 

Section El is one of the main sections of this paper. We first explain a general method on 
how one can efficiently design excellent constellations for any set of parameters M, N, T, L 
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and p. For this we review the properties of the complex Stiefel manifold and the Cayley 
transform. We conclude this section with an extensive table where we publish a large set 
of codes having some of the best diversity sums and diversity products in their parameter 
range. More extensive lists of codes with large diversity can be found on the website jTTIj . 

Finally in Section |H1 we explain how the algebraic structure which underlies most of the 
derived codes can be used to have a fast decoding algorithm. Our simulations indicate that 
in the design of codes more attention should be given to the diversity sum (more generally 
diversity function) which previously has not been fully studied. 

2 The Diversity Function, the Diversity Product (DP) 
and the Diversity Sum (DS) 

In this paper we will be concerned with the construction of constellations where the right 
hand sides in (jl.lj) and (jl-2j) . maximized over all pairs I, V is as small as possible for fixed 
numbers of T, M, N, L. As already mentioned this tasks depends on the signal to noise ratio 
the system is operating. For this purpose we define the exact diversity function dependent 
on the constellation V = {$1, . . . , $l} and a particular SNR p through: 

T> e (V,p) := maxProb( choose <Ev | transmitted ) (p) (2-1) 

For a particular constellation with a large number L of elements, with many transmit and 
receive antennas the function T> e (V,p) is very difficult to compute. Indeed for each pair 
it is required to compute the singular values of the M x M matrix and then 

one has to evaluate up to LiL — l)/2 integrals of the form (jl.lj) and this has to be done for 
each value of p. Although this task is formidable it can be done in cases where T, M, L are 
all in the single digits using e.g. Maple. 

Using Chernoff 's bound (jl.2j) we define a simplified function called the diversity function 
through: 

1 M 

P(V,p):=max-n 

m=l 

The computation of T>(V, p) does not require the evaluation of an integral and the compu- 
tation requires essentially the computation of ML(L — l)/2 singular values. The singular 
values 5 m (^&i') are by definition all real numbers in the interval [0, 1] as we assume that the 
columns of form both orthonormal frames. The functions T> e (V,p) and P(V, p) are 

the smallest if the singular values 5 m ($*<Ev) are as small as possible. These numbers are all 
equal to zero if and only if the column spaces of $/, <ly are pairwise perpendicular. We call 
such a constellation fully orthonormal. Since the columns of generate an M-dimensional 
subspace this can only happen if L < T/M. On the other hand if L < T/M it is easy to 
construct a constellation where the singular values of ($^$^) are all zero. Just pick LM 
different columns from a T x T unitary matrix. Figure Q depicts the functions V e (V,p) and 
T>(V, p) for a fully orthonormal constellation with T = 10 and M = N = 2. 
In order to study the function P(V, p) more carefully let 
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Figure 1: Diversity function T>(V, p) and exact diversity function T> e (V, p) of a fully orthonor- 
mal constellation. 

In some small interval [pi,p2] the maximum in (|2.2|) is achieved for some fixed indices 1,1' 
and in terms of p the function T>(V, p) is of the form: 



where the coefficients Cj, . . . ,Cm depend on the particular constellation and on the chosen 
interval [p\ , p^\ ■ For an interval close to zero the dominating term will be the coefficient c\ . Up 
to some factor this term will define the diversity sum of the constellation. When p » then 
the dominating term will be the coefficient cm and up to some scaling this term will define 
the diversity product of the constellation. A constellation will have a small diversity function 
for small values of p (and presumably performs well in this range) when the constellation 
is chosen having a large diversity sum. A constellation will have a small diversity function 
for large values of p (and presumably performs well in this range) when the constellation 
is chosen having a large diversity product. In the next two subsections we will study the 
limiting behavior of T>(V,p) as p goes to zero and to infinity. 
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2.1 Design criterion for high SNR 

When the SNR p is very large then T>(V, p) can be approximated via: 

ViV p) ~ max - ( ™ 2 TT l — 

P(V,p)_max 2 ^ 4(1 + pT/M) j 11 



(2.4) 



It is the design objective to construct a constellation $x, $ 2) • ■ • , $n such that 



rmnn(l-^($^0) 

is as large as possible. This last expression defines in essence the diversity product. In order 
to compare different dimensional constellations it is customary to use the definition: 

Definition 2.1. (See ^H]) The diversity product of a unitary constellation V is defined as 

i 

M 

2\ 



I]V = mm 

\m=l 



An important special case occurs when T = 2M. In this situation it is customary to 
represent all unitary matrices in the form: 

\/2 ( I 



**=-uj- (2 - 5) 

Note that by definition of the matrix f ^ is a M x M unitary matrix. The diversity 
product as defined in Definition 12.11 has then a nice form in terms of the unitary matrices. 
For this let A m be the mth eigenvalue of a matrix, then 

1 - = \\ m {2I M - *\*v - = \s 2 m (I M - = U 2 m {m v - ¥,). 

So we have 

m 1 m 1 

11(1- ^($*$,))w = _JJ 5m (^, _ = -| det(*,, - 

m=l m=l 

When T = 2M and the constellation V is defined as above, then the formula of the diversity 
product assumes the simple form: 

n v = \ <^< £ i det (^ - ( 2 - 6 ) 

We call a constellation V a fully diverse constellation if Yl V > 0. A lot of efforts have 
been taken to construct constellations with large diversity product. (See e.g. [TBI EH El HU 
123 EZl EH] ) • For the particular situation T = 2M with special form ()2.5|) the design asks for 
the construction of a discrete subset V = . . . , ^l} of the set ofMxM unitary matrices 
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U(M). When this discrete subset has the structure of a discrete subgroup of U(M) then the 
condition that V is fully diverse is equivalent to the condition that the identity matrix is the 
only element of V having an eigenvalue of 1. In other words the constellation V is required 
to operate fixed point free on the vector space C M . Using a classical classification result of 
fixed point free unitary representations by Zassenhaus [SI], Shokrollahi et al. [2BJ EZJ were 
able to study the complete list of fully diverse finite group constellations inside the unitary 
group U(M). Some of these constellations have the best known diversity product for given 
fixed parameters M, N, L. Unfortunately the possible configurations derived in this way is 
somehow limited. The constellations are also optimized for the diversity product and as 
we demonstrate in this paper for unitary space time modulation maybe attention should be 
given to the diversity sum. 

In most of the literature mentioned above researchers focus their attention to constel- 
lations having the special form ()2.5|) . Unitary differential modulation JTHJ is used to avoid 
sending the identity (upper part of every element in the constellation) redundantly. This 
increases the transmission rate by a factor of 2 to: 

_ log 2 (£) = 2 log 2 (£) 
M T 

Because of this reason we will also focus ourselves in the later part of the paper to the special 
form (|2.5jl as well. Nonetheless it will become obvious that the numerical techniques also 
work in the general situation. 

2.2 Design criterion for low SNR channel 

As we mentioned before a constellation with a large diversity sum will have a small diversity 
function at small values of the signal to noise ratio. This is particularly suitable when the 
system operates in a very noisy environment. When p is small, using Formula (J2.3|) . one has 
the following expansion: 

n [i + 4 n p + / pT/M) (i - = n [i + pa - *i0)] 

rn=l ^ III m=l 

M 

= i + p^(i-e(^^)) + o(p 2 ). 

m=l 

When p — > 0, i.e. p — > 0, we can omit the higher order terms 0(p 2 ) and the upper bound 
of , requires that 

^(1 - Sli^)) = (M - \\^ v f F ) 

m 

is large. In order to lower the pairwise error probability, it is the objective to make \\ F 
as small as possible for every pair of /, V . It follows that at high SNR, the probability primarily 
depends on [] t (l — <^($*<Ev)), but at low SNR, the probability primarily depends on 
£m=i In order to be able to compare the constellation of different dimensions, 
we define: 
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Definition 2.2. The diversity sum of a unitary constellation V is defined as 



5> 




M 



Again one has the important special case where T = 2M and the matrices take the 
special form ()2.5|) . In this case one verifies that 



12 i||/ + *r*HIJ 4 

1 , . , _ - - I 



-tr (2/ + + = -(4M - (2M - tr + %*%>))) 

= ~(4M - tr ((*, - - *,,))) = J(4M - ||*, - tt,,|lj0 

For the form (|2.5J) the diversity sum assumes the following simple form: 

VV = mm— =||tf,-tfp|L. (2.7) 
^ M' 2^M 

Without mentioning the term the concept of diversity sum was used in [T2j. Liang and 
Xia |2"T1 p. 2295] explicitly defined the diversity sum in the situation when T = 2M using 
equation (|2.7j) . Definition 12 .21 naturally generalizes the definition to arbitrary constellations. 

We want to point out that the diversity sum is the design criterion only for unitary 
constellation. Hochwald and Marzetta ^H] calculate the non-coherent space time channel 
capacity and indicate that unitary signal constellation are capacity achieving signal sets 
only for high SNR scenarios. For low SNR case, the transmitting power should be allocated 
unsymmetrically, i.e., unitary constellations are not capacity achieving in the first place. 
However unitary signal sets are easily manageable and one can take advantage of differential 
modulation technique [18 j to speed up the transmission. Moreover our simulation results 
indicate that codes with near optimal diversity sum tend to perform significantly better 
compared to the currently existing ones optimized for the diversity product for low and even 
moderate SNR scenarios. So it is quite reasonable and more toward the practical use to 
construct unitary constellations with good diversity sum. 

As the formulas make it clear the diversity sum and the diversity product are in general 
very different. There is however an exception. When T = 4, M = 2 and the constellation V 
is in the special (J2.5|) . If in addition all the 2x2 matrices {*!,..., *l} are a subset of the 
special unitary group 

SU{2) = {Ae C 2x2 | A* A = I and det A = 1} 

then it turns out that the diversity product F] V and the diversity sum ^ V of such a 

constellation are the same. For this note that elements of SU(2) have the special 
form: 

fa 6\ T ( c d 



Through a direct calculation one verifies that det(*i — = |a — c| 2 + |6 — <i| 2 and 
— ^z'IIf = 2(|a — c| 2 + \b — d\ 2 ). But this means that Y[V = Yl V f° r constellations 
inside SU{2). ^ 
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2.3 Four illustrative examples 

The diversity sum and the diversity product govern the diversity function at low SNR re- 
spectively at high SNR. Codes optimized at these extreme values of the SNR-axis do not 
necessarily perform well on the "other side of the spectrum". In this subsection we illus- 
trate the introduced concepts on four examples. All examples have about equal parameters, 
namely T = 4, M = 2 and the size L is 121 respectively 120. The first two examples are 
well studied examples from the literature. We derived the third and the fourth examples by 
geometrical design and numerical methods respectively. 

Orthogonal Design: This constellation has been considered by several authors j2j 128] . 
For our purpose we simply define this code as a subset of SU(2): 



The constellation has 121 elements and the diversity sum and the diversity product are both 
equal to 0.1992. 

Unitary Representation of SX 2 (F 5 ): Shokrollahi et al. [2H] derived a constellation using 
the theory of fixed point free representations whose diversity product is near optimal. This 
constellation appears as a unitary representation of the finite group SX 2 (F 5 ) and we will 
refer to this constellation as the SX 2 (F5)-constellation. The finite group SZ^Fs) has 120 
elements and this is also the size of the constellation. In order to describe the constellation 
let r\ = e~ and define 



Then the constellation is given by the set of matrices {PQ) J X, where j — 0, 1, • • ■ , 9, and X 
runs over the set 



{/ 2 , P, Q, QP, QPQ, QPQP, QPQ 2 , QPQPQ, QPQPQ 2 , 

QPQPQ 2 P, QPQPQ 2 PQ, QPQPQ 2 PQP}. 



The constellation has rate R = 3.45 and Y\SL 2 (¥ 5 ) = £SL 2 (F 5 ) = ±W ' ~ 0.3090. 



The diversity product of this constellation is truly outstanding. For illustrative purposes we 
plotted in Figure El the exact diversity functions and the diversity function of this constella- 



Numerically Derived Constellation: Using simulated annealing algorithm we found 
after short computation a constellation with very good diversity sum. The constellation is 
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Figure 2: Diversity function T>(V,p) and exact diversity function for the group constellation 
SL 2 (¥ 5 ). 



given through a set of 121 matrices 
tffe, :=A k B l \A = 



B 



-0.9049 + 0.3265 *i 0.1635 + 0.2188 * % 
0.0364 + 0.2707 *i -0.8748 + 0.4002 * i 

-0.1596 + 0.9767 * i -0.1038 + 0.0994 * i 
0.0833 - 0.1171 *i -0.9432 + 0.2995 *i 



,fc,Z = 0,1,- -- ,10 . 



As we explain in Section |Hl the maximum likelihood decoding of this constellation admits a 
simple decoding algorithm: sphere decoding. 



Geometrically Designed Constellation: Based on the algebraic structure we are going 
to propose in this paper, we further implement the geometrical symmetry into this structure. 
A geometrically designed constellation can be described as follows: 

/ „177r/60 f. 

% :=A k B k \A - 1 



Q e 137r/60 

/ cos(227r/60) sin(227r/60) \ 
V - sin(227r/60) cos(22tt/60) J ' 

This constellation has superb diversity sum and reasonably good diversity product. One can 
also use sphere decoding to implement maximum likelihood decoding of this constellation. 
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Table 1. The following table summarizes the parameters of the four constellations: 





Orthogonal 
design 


SL 2 (¥ 5 ) 


Numerically 
derived 


Geometrically 
designed 


Number of elements 


121 


120 


121 


120 


diversity sum 


0.1992 


0.309 


0.3886 


0.4156 


diversity product 


0.1992 


0.309 


0.0278 


0.1464 



Of course we were curious about the performances of these four different codes. Figure El 
provides simulation results for each of the four constellations. Note that the numerically 
designed code who has a very bad diversity product is performing very well nevertheless due 
to the exceptional diversity sum. One can see that up to 12db numerically derived codes 
outperform the group code by about 1 db. In fact, our simulation results show that until 
35db the numerical one is still performing much better than the orthogonal one. However at 
around 18db, the group constellation surpasses the numerical one due to exceptional diversity 
product. The geometrically designed constellation has better diversity sum and diversity 
product than the numerical one, therefore its performance is better than the numerical one 
(our results show that their performance curves are quite close, although the geometrical 
one is slightly better). These simulation results give an indication that the diversity sum is 
a very important parameter for a unitary constellation at low SNR regime. 




Figure 3: Simulations of four constellations having sizes T = 4, M = 2 and L = 120 
respectively L = 121. 
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3 Constellations With Algebraic Structure 



Before we venture into the realm of structured constellation, we would like explore random 
unitary space time constellations first. We introduce the Haar distributed random matrix, 
which in some sense can be viewed as a high dimensional generalization of a complex random 
variable with circular symmetric distribution CAf(0, 1). 

Definition 3.1. The Haar measure on U(M) is defined to be a probability measure Ti, on 
U(M) which is translate invariant: for any mensurable set S in U(M) and any fixed element 
U in U(M) 

H(S) =H(U S). 

A unitary random matrix U is Haar distributed (h.d.) if for any measurable set S we have 

Pr(U eS) = H(S). 

Remark 3.2. Note that h.d. matrix is also called isotropically distributed matrix in [22] • We 
want to point out that Haar measure can be defined more generally. In fact every compact 
Lie group admit a unique (up to scalar) translate invariant measure: Haar measure [3]. 

A well known yet non-trivial fact is that for any measurable set S C U(M), we have 

H(S)=H(S*), 

where S* consists of the conjugate transpose of all the elements in S. Thus for a h.d. matrix 
U, one can verify 

Pr(U* G S) = Pr(U E S*) = H(S*) = H(S). 

Immediately we conclude U* is also h.d. matrix. Also one can verify that the product of two 
h.d. matrices is still h.d. Another very interesting property about a h.d. matrix is about 
its spectrum. As derived in [Oj, the joint probability density for the eigenvalues of a h.d. 
random matrix U ~ diag (e^ 1 , e 102 , • • • ,e l9m ) in U(M) is given by the Weyl denominator 
formula: 

j<k 

The properties of h.d. matrices lead to the following theorem about random unitary space 
time constellation: 

Theorem 3.3. For a random unitary space time constellation V consisting of L h.d. inde- 
pendent random matrices Ui, U2, • • • , Ul, we have 

Pr(Y[V = 0) = 0, 

that is the probability ofV being fully diverse is 1. 
Proof. First we can rewrite 

Pr(JJV = 0) = Pr ( |J I det(Uj-Uk) | = 0) < Pr(| det(Uj - U k )| = 0). 

j<k j<k 
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Next we are going to show that the probability of the event | det(Uj — Uk)| = happening 
is 0. Now, 

Pr(| det(Uj - U k )| = 0) = Pr(| det(J - Uj*U k )| = 0). 

Let U denote Uj*Uk, we know U is h.d. matrix. Using the Weyl denominator formula, one 
computes 

M M „ „ 

Pr(| det(J - U) | = 0) = FrflJ ft = 0) < M £ / / II ^ ~ ^dfodfa ■ ■ ■ d9 M . 

i=i i ) ■ l=1 J J0i—0 j <k 



Since 



JI | e ^ - e idk \ de x de 2 ■ ■ ■ de M < 2 M ^ M ^ j j de 1 de 2 ■ ■ ■ de M = o, 

i< 

we conclude that 



e 3 =o j<k JJe 1= o 



Pr( | det(Uj -U k )| = 0) = 0. 



Consequently 



Pr(JJV = 0) = 0, 

that is the probability of V being fully diverse is 1. □ 

Note that if an M x M matrix G of independent complex Gaussian entries is input to 
QR algorithm, the resulting unitary matrix Q is Haar distributed [§]. For simplicity we 
sketch the proof as follows: First one can write Q = GR^ 1 , then for a fixed unitary matrix 
Uq, it can be checked that UqG has the same distribution as G. Consequently UoQ has 
the same distribution as Q, i.e., the distribution of Q is translate invariant. Therefore the 
uniqueness of translate invariant measure on a compact Lie group guarantees that Q is Haar 
distributed. As a consequence of the above theorem, an algorithm which produces a fully 
diverse unitary constellation with probability 1 can be given as follows: take L instance 
of complex Gaussian matrices and feed them through the QR algorithm, the resulting L 
unitary matrices constitute a fully diverse constellation with probability 1. 

From an algebraic geometry point of view one easily shows that the set of constellations 
with Y[ V = forms a lower dimensional proper algebraic sub- variety of U(M) L . In particular 
the set of all the fully diverse constellations is Zariski open [Tl| in U(M) L , i.e., fully diverse 
constellations are dense in U(M) L . Haar distributed random constellations won't be practical 
for maximum likelihood decoding in high transmission rate scenario because no algebraic 
structure is assumed for random constellation and therefore the decoding process will be too 
complex. In the sequel we are going to investigate structured constellations and explain how 
one can restrict the parameter space to judiciously chosen subsets and how one can convert 
maximum likelihood decoding to lattice decoding by using structured constellations. 

Consider a general constellation of square unitary matrices, 

V = {*i,tf 2 ,-- - 

In order to calculate the diversity product, one needs to do L ^ L 2 ~ 1 ^ calculations: | det(\l/j — 
for every different pair i, j. The same statement can be said about the diversity sum, however 
for simplicity we only show the diversity product case in the sequel unless specified otherwise. 



13 



If one deals with a group constellation then one needs only to calculate L — 1 such 
determinant calculations and this is one of the remarkable advantages of group constellations. 
This is a direct consequence of 

| det(*i - = I det(tfi) det(J - **^)| = | det(J - #*^)|, 

where is still in the group. 

As we mentioned before group constellations are however very restrictive about what the 
algebraic structure is concerned. In the following we are going to present some constellations 
which have some small number of generators and whose diversity can be efficiently computed. 
This will ensure that the total parameter space to be searched is limited as well. We start 
with an example: 

Example 3.4. Consider the constellation 

V = {A k B l \A,B G U(M),k = 0, • • • ,p,Z = (),••• ,q}. 

The parameter space for this constellation is U(M) x U(M), this is a manifold of dimen- 
sion 2M 2 and the number of elements in V is (p + l)(q +1). If one has to compute 

|det(\I / j — for every distinct pair this would require ( ^ ^ ^ ^ determinant 

calculations. We will show in the following that the same result can be obtained by doing 
2pq + p + q determinant computations. 

Let and be two distinct elements having the form A kl B lx and A h2 B l2 respectively. 
We have now several cases. When k\ = k 2 , then necessarily li ^ l 2 and the distance is 
computed as 

| det(A kl B h - A k2 B h )\ = | det(J - B\ h ~ ll % 

where | Z 2 — Z x | is an integer between 1 and q. If li = l 2 , then we have k\ ^ k 2 and the distance 
is computed as 

| det^ 1 ^ 1 - A k2 B h )\ = | det(J - A lk2 ~ kll )\, 

where \k 2 — ki\ is an integer between 1 and p. If {ki < k 2 and l\ < l 2 ) or {ki > k 2 and l\ > 
l 2 ), we have 

| det(A kl B h - A k2 B l2 ) \ = | det(J - A^ k2 ~ k ^ B^ h ~ 1 ^) , 

where 1 < \k 2 — k x \ < p and 1 < \l 2 — h\ < q. Similarly if (/q < k 2 and / x > l 2 ) or 
{k\ > k 2 and Z x < l 2 ) then 

| det^ 1 ^ 1 - A k2 B l2 )\ = | det(A |fe2 - fel1 - B^ l ^)\, 

with 1 < \k 2 — ki\ < p and 1 < | — Zi j < p. The total number of distances to be computed 
is in total equal to 2pq + p + q. 

The number of distances to be computed indicates how complex the calculation for the 
diversity is. In fact the smaller this number is, intuitively the larger possibility of finding 
a unitary constellation with good diversity we will have. An immediate observation is that 
for two pair of unitary matrices (A, B), (C,D), if (C, D) = (UAV,UBV) or (C, D) = 
(UA^ViUB^V), then | det(C — D)\ = | det(A — B)\. We are going to consider several 
constellations starting with this observation. 
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Example 3.5. Consider the case that G C U (n) is a subgroup with L elements, then for any 
two distinct elements A, B e G we have | det(A - B)\ = \ det(J - A~ X B)\ with A~ l B E G. 
Therefore at most L — 1 distance calculations are needed to derive the diversity product. 
The product of two group constellations has the similar property. Consider Gi C\U (M) with 
order l iy where i = 1,2. Let 

G = {AB\A eG u Be G 2 }. 

Since | det(A XJ Bi - A 2 B 2 )\ = | det(I - A^A 2 B 2 B^)\ with A^A 2 B 2 B^ E G, at most L - 1 
calculations are needed in this case, where L = l\l 2 . 

Example 3.6. Consider a constellation with the following form: 

{j£B*\i = 0,-- • ,h - l,j = 0,- •• ,Z 2 - 1 and A, Be U(M),A h = I,B h = I}. 

It can be checked that for the above constellation at most L — 1 calculations are needed, 
where L = l\l 2 . 

Group structures do have certain advantages for constructing unitary constellations: it 
is less complex to calculate the diversity product (or sum); the possibility of finding a large 
diversity constellation intuitively may be increased. However the constellations found by this 
approach [2E] are really few and far between. Somehow one wonders if the group structure 
is too restrictive to find a good-performing constellation. 

In the sequel we are going to loosen the constraints imposed by the group structures. 
As demonstrated in Example 13 . 41 it is desirable to have a small dimensional manifold (in 
Example 13.41 it was U(M) x U(M)) which parameterizes a set of potentially interesting 
constellations. Having such a parameterization will help to avoid the problem of "dimension 
explosion". The set of constellations parameterized by U(M) x U(M) in Example 13.41 are 
interesting as we are not required to compute all pairwise distances in order to compute the 
diversity product (sum). 

Definition 3.7. Let X be the set {xi,x 2 , • • • , x n } and F be the free group on the set X. A 
subset G C U(M) is called freely generated if there are elements {gi,g 2 , ■ ■ ■ ,g n } C G such 
that the homomorphism : F — > G with <f>{xj) = g^ is an isomorphism. 

An immediate consequence of this definition is that every element in G can be uniquely 
written as a product of <?j's and g^ lj s. The elements are called the generators of G. A 
freely generated subset G is simply parameterized by the set: 

{a^af 2 • • • a^ k | a$ is one of g-s, pi e Z} . 

Take an element g G G with its representation g = Yi i=1 af% we say that the presentation 
is reduced whenever a« ^ a i+ i for i = 1, ... ,n — 1. Observe that taking the product of 
distinct matrices Yl7=i ^ ^ s numerically expensive, however taking the power of one matrix 
A k is much easier (note that for A = U U^ 1 with diagonal, we have A k = U U^ 1 ). 
Moreover by considering the powers of one matrices, we are able to impose the lattice struc- 
ture to the constellation, which makes sphere decoding of structured constellations possible, 
(see Section IHI) Therefore we are interested in "normal" elements of G. 
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Definition 3.8. We say that an element g = Yli=i a T in reduced form is a normal element 
whenever a, 7^ aj for i 7^ j. A subset V of the freely generated set G is said to be a normal 
constellation if every non-identity element in V is normal. 

Since finding an inverse of a matrix is numerically expensive, we also limit our searches 
to positive constellations: 

Definition 3.9. An element g in G with the reduced form g = Y\i=i a i Pl i s sa id to be a 
positive element if Pi > for i — 1, 2, • • • , k. A subset V of the freely generated set G is said 
to be a positive constellation if every non-identity element in V is positive. 

Positive normal constellations are desirable for numerical searches as they can be effi- 
ciently parameterized and searched. If one wants to compute the diversity product (or sum) 
of an arbitrary positive constellation with L elements one still has to compare a total of (J) 
pairs of matrices. In the sequel we will impose more structure on a constellation VcG which 
will guarantee that only L — 1 pair of elements have to be compared during the diversity 
product (sum) computation. 

Definition 3.10. Two unitary matrices A, B G G are said to be equivalent (denote by 
A ~ B) if there is unitary matrix U G G such that A = UBU' 1 or A = UB^U^ 1 . [A] 
will denote all the matrices that are equivalent to A. For a constellation V C G, we say 
V = {^1, ^2, ' ' ' , ^l} has a weak group structure if for any two distinct elements the 
product is equivalent to some ^fc. 

The reader verifies that we indeed defined an equivalence relation. Note also that V has 
a group structure as soon as ^^^i ^ s always another element of V and this explains our 
wording. 

Lemma 3.11. Let V = {^o — 1,^1,^2, '• • be a constellation with a weak group 

structure. In order to compute the diversity product (sum) it is enough to do L — 1 distance 
computations. 

Proof. 

I det(^ - = I det(J - ttr 1 **)! = I det ( J - B )l 

where B e V is an element in V equivalent to This shows the result for the diversity 

product. If one is concerned with the diversity sum then the same argument still holds if the 
absolute value of the determinant | det(-) | is replaced by the Frobenius norm || . \\ F . □ 

Based on this lemma we are interested in finite constellations inside G whose elements 
have a weak group structure and are all normal. The following theorem provides a complete 
characterization of all these constellations: 

Theorem 3.12. Let V C G be a finite positive normal constellation (including identity 
element) with L > 3 elements. If V has a weak group structure then V takes one of the 
following forms: 

. {I,A,A\...,A L ~i} 

• {I, AB, A 2 B 2 , ■ ■ ■ , A L ^ 1 B L ^ 1 } 
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where A = gf' 1 , B = g Pj for some i ^ j. 

The proof of Theorem 13.121 is rather involved. In order to make it more understandable 
we will divide it in several definitions and lemmas. 

Definition 3.13. For any element $ G G, we define the length of \I/ = Yli=i a i Pi t° be 

k 

length (\&) = / ]pi- 

i=l 

It is a routine to check that the definition is well-defined and doesn't depend on the 
representation of the element. For the identity element one will have length (/) = 0. One 
immediate consequence from this definition is that if A ~ B, one will have |length(A)| = 
|length (B)\. The following lemma claims that any freely generated positive weak group 
constellation "approximately" takes cyclic form. 

Lemma 3.14. Let V = {^>q = I, ^i, \&2, • ■ ■ > ^l-i} C G be a positive constellation of the 
freely generated set G C U(M). Suppose length < length (tyj) for i < j. If V is a weak 
group constellation, then 

where = {aia 2 ■ ■ ■ Oi|ai, a 2 , • • • , <2j G [^i]}- 

Proof. We first show that length (^) < length i^j) for i < j: Indeed, if length = 
length (\F,), then length (v^ 1 ^) = length - length (^) = 0. That means ^^^j ~ /, 
equivalently one will have ty^iffj = /, i.e. = tyj. That contradict the fact that and 
tyj are distinct. 

Consider ^^ 1 ^ 2 - Since < length (Nt^ 1 ^) = length (* 2 ) - length (^i) < length (9 2 ), 
therefore ^T 4 ^ = ^i where ~ ^f x . So \I> 2 = ^i^i G [^l] 2 - Proceed by induction, one 
can show ^ 1 ^ k+1 = # 2 where # 2 ~ ^i. So = ^ fc ^ 2 G by induction. □ 

Remark 3.15. An immediate observation is that 

length = i * length (*i). 
Take two positive normal elements in G with their reduced forms: 

We define the shift operator S k on the reduced form of a positive normal element \I/ by 
induction: S^) = Si (af a 2 2 ■ ■ ■ = a p 2 2 ■ ■ ■ a^pa^ 1 and S k +i = S k o Si. We assume that 
5o( l I / ) = then apparently for a fixed element \I/ shift operator is periodic. We have the 
following lemma. 

Lemma 3.16. ~ \I/ 2 if and only if = S k (^2) for some k. 
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Proof. The sufficiency part of this lemma is straightforward. So we have to prove the ne- 
cessity part. Since ~ ^2, according to the definition of equivalence there exists c such 
that c^ic~ l = ^2 or c^ic -1 = ^2 ■ However since length (c^ic -1 ) = length (^2) > and 
length (^Z^ 1 ) < 0) the second case won't happen. The only possibility is c^icT 1 = ^2- We 
assume that c is generated by only one generator and further assume c = c^ 1 with l\ > 0, 
then we will have 

c l lafaf---a^ h = bfbf.--b%. 

So Ci = a m and l\ < p m follows, otherwise the left hand side of the equation above will 
have negative power, while the right hand side only has positive power. This will contradict 
the uniqueness of the representation of the same element. In fact li = p m , since otherwise 
^2 — c^a^a^ 2 ■ ■ •c^ m_/l . This will contradict the fact that ^2 is a normal element. So with 

n P m n Pl ■ ■ ■ (7 Pm_1 — h qi h q2 ■ ■ ■ h q " 
u m u l a m-l — u l u 2 u n ; 

one can check m = n and ^2 = 5'm-i(^ / i)- 

Proceed by induction, suppose c has the reduced form c = c^Cg • • • c^tii then the follow- 
ing equation follows: 

Jl J2 J k+1 n Pl n P2 n Pm„- l k+l „-h„-h _ 7,917,92 Wn 
c l c 2 c k+l a l a 2 a m L k+1 C 2 c l — °1 u 2 u n • 

Without loss of generality, we assume lk+i > and apply the same argument as in the one 
generator case. One proves a m = Ck+i and 4+1 = p m . Therefore we reach the following 
equation: 

c'/cff ■ ■ ■ c^S^^Jc- 1 * ■ ■ ■ c 2 h c- h = bfbf 
By induction, \l>2 = Sfci S'm-it^i) = 

+ m _i(\l/i) for some k%. O 

Proof of Theorem \S.lSl Pick any two distinct elements G V having length (^i) < 

length (tyj). We claim that if = a\CL2 • • • a m , then either there exists 1 < k < m — 
1 such that ^/j = a\a2 • • • akbib 2 ■ ■ ■ hdk+i ■ ■ ■ a m , or = 6i6 2 • • • &«Oia 2 • • • a m or = 
aia 2 • • • a m bib2 ■ ■ - h for some / > 0. 

Suppose that the claim is not true, then for \^ = C1C2, • • • c p , there exist k\, k 2 such that 
< ki < m, 1 < k 2 < m + I and k\ < k 2 — 1 and \&j will take the following form: 

vj/j- = ai a 2 ■ ■ ■ a kl bib2 ■ ■ ■ ka k2 ■ ■ ■ a m , 

where b x 7^ a kl+ i and b\ 7^ a k2 -\. (For the special case k x = 0, we assume C\ 7^ a\. For 
the special case k 2 = m + 1, we assume c p 7^ a m .) Then would be equivalent 

to aZ_ x • • ■ dkl + ib\b2 • • - bi, which in any case won't be equivalent to any positive element 
= did2 ■ ■ ■ d q or /. That contradicts the fact that V is equipped with a weak group 
structure. 

As explained above we can further assume that 

length (/) < length < • ■ ■ < length (^l-i). 

If is generated by only one generator, i.e. ^1 = gf' for some i. Since ^2 is a normal 
element, according to the claim, either ^2 = ^1^2 or ^2 = ^2^1 for some ^2- In either 
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case ^2 will be equivalent to ^i, while Lemma 13.161 will guarantee ^2 = 3>i- Therefore we 
will have \l/ 2 = 9i Pi ■ Proceed by induction, it can be checked that ^1 = gj Pi for every /. So 
the constellation will take the first form in the theorem. 

If $1 is generated by two generators, i.e. \&i = g Pt g P3 for some i,j. According to the 
claim, we will have ^2 = ^1^2 or ^2 = ^2^1 or ^2 = <?f l ^2g/ ■ Because ^2 is equivalent to 
^1, ^2 is a shifted version of ^1. Exhausting all the possibilities, the first two cases would 
make ^2 a non-normal element, so the only possibility is the third case. Consider two shifted 
version of ^j: Sq^i) = 9^9^ an d = g Pi ' gf* . Only 5o(^i) will satisfy the condition 

that ^2 is a normal element. So the analysis above shows that 

*2 = g?*ifl = g? l gT- 

By induction it can shown that 

v/c+i - g\ Vkgj - gi gj 

So in this case, the constellation will take the second form in the theorem. 

However the constellation doesn't exist if is generated by more than 3 elements. 
Indeed suppose with the reduced form ^ = a^a?, 2 ■ ■ ■ aP^ with m > 3, then ^ 2 will take one 
of the following form: \j> 2 af af " " " a £T> a f ^2^2 " " " aP mi " " " > a? " " " <T^2 with ^ 2 being 
a shifted version of But ^2 wouldn't be a normal element for any of the above form, so 
there doesn't exist weak group constellation for this case. 

□ 

A weak group constellation is very group like, while it is not exactly a group. It does 
keep the advantage of a group constellation: for example, for any weak group constellation 

V taking the second form in the theorem, only L — 1 computations | det(J — A k B k )\ for 
k = 1,2, ■ ■ ■ ,L — 1 are needed to calculate the diversity product. It also overcome the 
disadvantage of group codes: one can freely choose the generators, while in group structures, 
the generators have to satisfy certain relations to be a group. Last but not least it turns 
out that the restriction to code elements in normal form is very advantageous during sphere 
decoding. In the next section we will mainly use the second weak group structure as described 
in Theorem 13. 121 Before we describe these search procedures we would like to illustrate some 
alternative methods. 

It is possible to increase the number of generators to obtain new structures. For instance, 

V = {A k B l C m \A,B,C EU(M),k = 0,--- ,p, I = 0, • • • , q, m = 0, • • • , r}. 

For a unitary constellation V = = !,■■■, L}, we call V s = {U§iV\i = 1, • • • , L} 

shifted version of V. It will be straightforward to prove that V s has the same complexity 
as V when one calculates the diversity. {A k CB k \A, B,C G U(M),k — 0, • • • , L — 1} is a 
shifted copies of the second weak group structure in Theorem 13.121 To see this, note that 
A k CB k = A k CB k C- 1 C = A k (CBC- 1 ) k C. It can checked that A k B L+1 ~ k = A k (B- 1 ) k B L+ \ 
therefore {A k B L+1 ~ k \A, B G U(M),k = 1, ■ ■ ■ , L} is also a shifted version of the second 
form weak group structure. 

Also we can consider the "combination" or the "product" of two structures. For ex- 
ample, {J, A, AB, ABA, ABAB, ABABA, ■■■} is the union of {(AB) k \k = 0, • • • } and its 
shifted version {(AB) k A\k = 0,---}. Another example is the product case: let Vi = 
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{/, C, C 2 , C 3 , • • • } and V 2 = {I, A, AB, ABA, • • • } and consider the Cartesian product con- 
stellation 

V = Vi x V 2 = {AB\A eVi,Be V 2 }. 

One may wonder how restrictive the proposed structures are. We all know a compact 
Lie group can be generated by any open neighborhood of any element in the Lie group. 
So with the above structure, even if one chooses the generators locally, the elements in the 
constellation could be spreading out on the whole manifold. Somehow this indicates that 
the proposed structure won't be too restrictive. 

4 Geometrical Design of Unitary Constellations with 
Good Diversity 

For low dimensional constellations, one may further specify the generators in the proposed 
structure. Observe that for the second form weak group constellation, one can always assume 
A is diagonal. In the sequel, we further assume that B is real orthogonal, i.e. based on the 
weak group structure we consider the following 2 dimensional constellation: 

V = {A k B k \A=( e ™ l),B=( C ° SZ SinZ ),k = 0,l,..- ,L-1}. (4.1) 
1 1 V e w J V — sin z cos z J J v ' 

There are several ways to design constellations with good diversity from this specific 
structure. A natural idea is to do Brute Force search using fine step size. Another approach 
is to design the constellation with the help of geometrical intuition. Note that a 2 x 2 com- 
plex matrix can be viewed as a vector in C 4 . In this context A and B can be viewed as 
"rotation" transforms (induced by regular matrix multiplication) acting on C 4 . A constel- 
lation of form (|4.1j) can be viewed as a set of rotated vectors under the transforms A h B k , 
k = 0, 1, • • • , L — 1. Intuition says that good constellations can be found if the rotation 
angle is symmetrical. Based on the idea above we assume that x, y, z to be the multiples of 
2%/L, we found a lot of good codes resulted from this geometrical symmetry (see tables in 
Section IHJ). 

2 dimensional constellation design has been studied in |2T] . In this paper Liang proposed 
very interesting parametric codes and many codes with excellent diversity are found. The 
codes shown in [2T] can be achieved by our design as well. In fact, most of Liang's codes 
belong to a special form of our parameterization (j4.1|) . To our best knowledge, most of our 
codes shown on the web site ^H] are the best codes ever found or never found before. 

Example 4.1. A very interesting code with 120 elements is found using this approach: 

v _ fAkf>k lA _ ( e n/m \ / cosvr/4 sinvr/4 \ ,_ ni nql 

V-{AB\A-^ Q elWm J,5-(^_ sin7r/4 cos?r/4 J,*-0,l,...,119}. 



It can be checked that Y\V = = \\J — 2> ^ ne diversity product and the 

diversity sum are identical to the ones of the S , L 2 (F 5 )-constellation. We simulated the per- 
formance of this code and compared it with the performance of the SX 2 (F 5 )-constellation. To 
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Comparisons of Two Constellations 




-8 -6 -4 -2 2 4 6 8 10 

SNR(db) 



Figure 4: 2 dimensional weak group constellations and group constellation 

our big surprise our new code performed considerably better than the 5X2(F5)-constellation. 
The constellation V with sphere decoding outperformed the SX2(F 5 )-constellation by about 
ldb up to about 20db (see Figure HJ. As the SNR goes higher, the two curves are coming 
closer though. 

In order to understand the difference in the performance of the two seemingly similar 
constellations we investigated the diversity product (DP) and diversity sum (DS) distance 
spectrum for each of them. As we explained before, for a unitary constellation with L 
elements, L(L — l)/2 distance calculations may produce distances with multiplicities. For 
example consider V as above, 360 out of 7140 pairs of elements have distance 0.3090 (see DP 
distance spectrum in Table EJ). So one can explain the behavior difference of the two codes 
using their distance spectrum. The following Table El shows that the DP and DS distance 
spectrum of our weak group constellation. 
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Table 2. 



Weak group constellation 
DP distance spectrum 



Weak group constellation 
DS distance spectrum 



distance 


distribution 


0.3090 


360 


0.3136 


480 


0.3895 


480 


0.3931 


1440 


0.4402 


240 


0.5000 


120 


0.5878 


120 


0.6360 


1440 


0.6787 


480 


0.7071 


600 


0.8090 


360 


0.8430 


480 


0.8660 


120 


0.8979 


240 


0.9511 


120 


1 


60 



1 • j 

distance 


distribution 


0.3090 


120 


r\ a a r\ c~\ 

0.4402 


240 


0.5000 


120 


0.5023 


480 


0.5457 


240 


0.5878 


120 


0.6367 


480 


0.6502 


240 


0.7071 


3000 


0.7598 


240 


0.7711 


240 


0.8090 


120 


0.8380 


240 


0.8647 


480 


0.8660 


120 


0.8979 


240 


0.9511 


120 


1 


60 



One can check that the DP distance spectrum of the SX2(F5)-constellation is identical to 
the DS distance spectrum. The following Table El shows that the DS distance spectrum for 
the S , L 2 (F 5 )-constellation has denser small distance distribution compared to DS spectrum 
of our constellation and this explains the considerable worse performance of this constellation 
in our simulations. 



Table 3. 



SL 2 (F 5 )-constellation 
DP (DS) distance spectrum 



distance 


distribution 


0.3090 


720 


0.5000 


1200 


0.5878 


720 


0.7071 


1800 


0.8090 


720 


0.8660 


1200 


0.9511 


720 


1 


60 



Although we have concentrated so far in the design of 2-dimensional constellations there 
is actually no restriction with our approach. The similar "rotation" idea can be applied to 
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other low dimensional constellation design. For instance, we can make further specifications 
to 3 dimensional weak group constellations: 



cosx sinx \ / e lz 

V = {A k B k \A=\ -sinx cosx ,5= cosw sinw | , k = 0, 1, • • • , L-l}. 

e iy J V — sin w cos w 

where x, y, z, w is assumed to take the multiple of 2n/L. Apparently algebraic design based 
on geometrical symmetry can be applied to any other structure as well. For instance consider 
the following specified structures: 

V = {A k B l \A=( 6 ™ l),B=( C ° SZ SinZ ),k = 0,l,...,p-l,l = 0,l,...,q-l}. 

L 1 V e y J V — sin z cos z J y 

where we can take x, y to be multiple of 2ft /p and z to be multiple of 2ir/q. We refer to 10 
for the designed low dimensional constellations from these approaches. 



5 Numerical Design of Unitary Constellation with Good 
Diversity 

In order to numerically design constellations, it will be necessary to have a good param- 
eterization for the set of unitary constellations having size L, operating with M transmit 
antennas. In this section we show how one can use the theory of complex Stiefel manifolds 
and the classical Cayley transform to obtain such a parameterization. 

5.1 The complex Stiefel manifold 
Definition 5.1. The subset of T x M complex matrices 

S TM := {$ G C TxM | $*$ = I M ) 

is called the complex Stiefel manifold. 

From an abstract point of view a constellation V := {$i, . . . , $l} having size L, block 
length T and operating with M antennas can be viewed as a point in the complex manifold 

M. := (St,m) L = $t,m x • • • x St,m ■ 

■ v ' 

L copies 

The search for good constellations V requires hence the search for points in M. whose diversity 
is excellent in some interval [p\,p2\- 

Stiefel manifolds have been intensely studied in the mathematics literature since their 
introduction by Eduard Stiefel some 50 years ago. A classical paper on complex Stiefel 
manifolds is [3], a paper with a point of view toward numerical algorithms is [7]. The major 
properties are summarized by the following theorem: 
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Theorem 5.2. St,m is a smooth, real and compact sub-manifold of C = M. of real 
dimension 2TM — M 2 . 

Some of the stated properties will follow from our further development. The following 
two examples give some special cases. 

Example 5.3. 



is isomorphic to the 2T — 1 dimensional unit sphere S . 

Example 5.4. When T = M then St,m = U(M), the group of M x M unitary matrices. It 
is well known that the Lie algebra of U(M), i.e. the tangent space at the identity element, 
consists of all M x M skew-Hermitian matrices. This linear vector space has real dimension 
M 2 , in particular the dimension of U(M) is M 2 as well. 

A direct consequence of Theorem 15.21 is: 

Corollary 5.5. The manifold M. which parameterizes the set of all constellations V having 
size L, block length T and operating with M antennas forms a a real compact manifold of 
dimension 2LTM - LM 2 . 

As this corollary makes it clear a full search over the total parameter space is only possible 
for very moderate sizes of M, L, T. It is also required to have a good parameterization of 
the complex Stiefel manifold St,m and we will go after this task next. 

The unitary group is closely related to the complex Stiefel manifold and the problem of 
parameterization ultimately boils down to the parameterization of unitary matrices. For this 
assume that $ is a T x M matrix representing an element of the complex Stiefel manifold 
$t,m- Using Gramm-Schmidt one constructs a T x (T — M) matrix V such that the T x T 
matrix [$ | V] is unitary. Define two T x T unitary matrices [$i | Vi] and [$2 I Vz] to be 
equivalent whenever $1 = $ 2 - A direct calculation shows that two matrices are equivalent 
if and only if there is (T — M) x (T — M) matrix Q such that: 



Identifying the set of matrices Q appearing in ()5.1|) with the unitary group U(T — M) we 
get the result: 

Lemma 5.6. The complex Stiefel manifold St,m is isomorphic to the quotient group 





(5.1) 



U(T)/U(T — M). 



This lemma let us verify the dimension formula for St,m stated in Theorem 15.21 



dimiS 



'T,M — 



dim U(T) - dim U(T — M) = T 2 — (T — M) 2 = 2TM - M 2 . 



24 



The section makes it clear that a good parameterization of the set of constellations 
V requires a good parameterization of the manifold M. and this in turn requires a good 
parameterization of the unitary group U(M). 

Once one has a nice parameterization of the unitary group U(M) then Lemma pro- 
vides a way to parameterize the Stiefel manifold St.m as well. Parameterizing U(T) modulo 
U(T — M) is however an 'over parameterization'. Edelman, Arias and Smith [7j explained a 
way on how to describe a local neighborhood of a (real) Stiefel manifold St,m- The method 
can equally well be applied in the complex case. We do not pursue this parameterization in 
this paper and leave this for future work. 

In the remainder of this paper we will concentrate on constellations having the special 
form (|2.5|) . From a numerical point of view we require for this a good parameterization of 
the unitary group and the next subsection provides an elegant way to do this. 

5.2 Cayley transformation 

There are several ways to represent a unitary matrix in a very explicit way. One elegant way 
makes use of the classical Cayley transformation. In order that the paper is self contained 
we provide a short summary. More details are given in [23 Section 22] and [To] . 

Definition 5.7. For a complex M x M matrix Y which has no eigenvalues at —1, the Cayley 
transform of Y is defined to be 

Y e = (I + Y)-\l-Y), 
where I is the M x M identity matrix. 

Note that (/ + Y) is nonsingular whenever Y has no eigenvalue at -1. One immediately 
verifies that (Y c ) c = Y. This is in analogy to the fact that the linear fractional transformation 
f( z ) = irf nas the property that f(f(z)) = z. Recall that a matrix M is skew-Hermitian 
whenever A* = —A. The set ofMxM skew-Hermitian matrices forms a linear subspace 
of C MxM = M. 2M having real dimension M 2 . This is the Lie algebra of the unitary group 
U(M). The main property of the Cayley transformation is summarized in the following 
theorem. (See e.g. [E, 26 ). 

Theorem 5.8. When A is a skew-Hermitian matrix then (I + A) is nonsingular and the 
Cayley transform V := A c is a unitary matrix. Vice versa when V is a unitary matrix which 
has no eigenvalues at —1 then the Cayley transform V c is skew-Hermitian. 

This theorem allows one to parameterize the open set of U(M) consisting of all unitary 
matrices whose eigenvalues do not include —1 through the linear vector space of skew- 
Hermitian matrices. The Cayley transformation is very important for the numerical design 
of constellations because it makes the local topology of U(M) clear. One can see that most 
optimization method require us to consider the neighborhood of one element in U(M). 

5.3 Simulated Annealing (SA) Algorithm 

In our numerical experiments we have considered several methods. Because there are a 
large number of target functions the best known optimization algorithms such as Newton's 
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Methods j2H H] and the Conjugate Gradient Method j2H U\ are difficult to implement. 
Surprisingly the Simulated Annealing Algorithm turned out to be very practical for this 
problem. 

Simulated Annealing (SA) is a method which mimics the process of melted metal getting 
cooled off. In the annealing process of the melted metal, first the metal is heated to melt, 
then the temperature is getting down gradually. The metal will get to a minimized energy 
state if the temperature is lowering slow enough. For more details about this algorithm, we 
refer to [H EB E5| ■ 

In fact, we would rather call it a general method instead of a concrete algorithm. Gen- 
erally speaking, for a given optimization problem we always take an initial solution in some 
certain way, then consider a second solution in the "neighborhood" of this solution. We will 
accept the solution according to some predefined criterion which might involve a probability 
threshold. 

Combining with good algebraic structure and Cayley transform, which is a good rep- 
resentation of any dimensional unitary matrix, one can see that numerical method can be 
applied to any dimensional and any size constellation design. Our implementation of the 
algorithm can be summarized in the following way, one can find simple sample program on 
our web site [TO] . 

1. Choose a proposed algebraic structure for the constellation. 

2. Generate initial generators of the whole constellation. One can either take an existing 
constellation as the start point or just take the initial point randomly. 

3. Generate randomly a new constellation using Cayley transform in the neighborhood 
of the old constellation where the selection is done using a Gaussian distribution with 
decreasing variances as the algorithm progresses. 

4. Calculate the diversity function (product, sum) of the newly constructed constellation. 

5. If the new constellation has better diversity function (product, sum), then accept the 
new constellation. If not, reject the new constellation and keep the old constellation 
(or accept it according to Metropolis's criterion [23]). 

6. Check the stopping criterion, if satisfied, then stop, otherwise go to 2 and continue the 



Example 5.9. As we mentioned before, one can either choose an existing constellation as 
the starting point for our numerical method or just take the initial point randomly. In the 
sequel, we use the group constellation (^21,4 in (2SJ: 



iteration. 




One can verify that 



JJVi = 0.3851 
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It seems that £21,4 is already a very good constellation, our algorithm only improves a little 
(see V2 below). However one can check for most of the cases, the algorithm will improve 
much compared to the original group constellation. 



where 



V 2 = {A k B l \k = 0, 1, • • • , 20, 1 = 0, 1, 2}, 



0.9415 + 0.3155 *i 0.0573 - 0.0222*2 0.0496 + 0.0882*? 

A=| 0.0160 - 0.0555 *i 0.4005 + 0.9136 * % 0.0326 - 0.0212 * % \. 

0.0579 + 0.0855 *i -0.0312 - 0.0099 * i 0.1384 - 0.9844 * i 

0.0175 + 0.0095 *i 0.9997 + 0.0111 *2 0.0079 + 0.0042*2 

D=\ 0.0086 + 0.0100 *i -0.0082 + 0.0040 * % 0.9999 + 0.0036 * % 

-0.4836 + 0.8750*? 0.0004 - 0.0198 * i -0.0045 - 0.0126 * i 



One verifies that 



Y[v 2 = 0.3874. 



Example 5.10. Different industrial applications require different level of reliability of the 
communication channels. One may want to optimize the constellation at certain Block 
Error Rate (BER) or Signal Noise Ratio (SNR). It can also be shown theoretically that 
numerical methods together with the proposed structure works in the same way if one wants 
to optimize the diversity function at a certain SNR. This is essential the case because for a 
complex matrix A and unitary matrices U, V one has that 

6 m (UAV)=6 m (A), (5.2) 

for m — 1, 2, • • • , M. With the constellation structures as above we are able to reduce the 
dimension of the parameter space and at the same time we have a considerable reduction 
in the number of targets to be checked. Intuitively algebraically designing codes for this 
purpose seems to be impossible. 

The following graph shows the comparison of three constellations with different dimen- 
sions with 2 receiver antennas. The first one is a 2 dimensional constellation with 3 elements 
(R = 0.7925) and optimal diversity product 0.8660 and optimal diversity sum 0.8660. The 
second constellation is a 3 dimensional constellation which has 5 elements (R = 0.7740) with 
diversity product 0.7183 and diversity sum 0.7454. The third constellation is a 4 dimensional 
one consisting of 9 elements (R = 0.7925) with diversity product 0.5904 and diversity sum 
0.6403. Here based on the structure A k B k we used Simulated Annealing to optimize the 
diversity function at 6db to acquire the last two constellations. 

One can see that around 5 db, the second constellation surpasses the first one and is 
getting better and better as the SNR becomes larger. This can be easily understood since 
the diversity function of the first constellation is approximately dominated by 1/p 4 at high 
SNR, while the diversity function of the second constellation is dominated by 1/p 6 . The 
same explanation can be applied to the third constellation's performance. One can even 
foresee that higher dimensional constellations will perform even better and the BER curve 
will be sharper than the lower dimensional ones. It is believable that higher dimensional 
constellations will achieve much more diversity gain compared to lower dimensional ones. 
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Comparison of Different Dimensional Constellation 
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Figure 5: Performance of different dimensional constellations with the same rate 



Surprisingly SA works very well when it is applied to an algebraic structure with symme- 
try Like all the other numerical methods, one has to suffer the loss of performance due to the 
increasing complexity as the size and dimension go up and due to the limited computational 
resources. However without any doubts, the numerical approach is very flexible and can be 
used for any dimensional and any size constellation, producing very good diversity. So a lot 
of good-performing unitary constellations are found this way, which were never found by any 
algebraic method. At the end of this section we will show some 2 dimensional constellations 
we found using various methods based on the proposed structure. We skip, however, our 
numerical results on the higher dimensional unitary constellation design, since one can check 
them on the web site [TU] . 

One very interesting fact is the numerical results for diversity sum from 3 dimensional 
structured constellation are even better than the corresponding upper bound for 2 dimen- 
sional constellations. Somehow it won't be too surprising if one notices that from U(2) to 
U(3), we have 5 more dimensions to manoeuvre. 

In [T3] packing problems on compact Lie groups are analyzed and the upper bound for 
the diversity sum and the diversity product are derived. In the following figure one can 
see the limiting behavior of 2 dimensional structured constellations compared to the upper 
bound. One can check for the comparisons for other dimensions. 
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Figure 6: 2 dimensional weak group constellations and upper bound 
5.4 Constellations with extremely large diversity 

In this subsection we list the best 2- dimensional constellations we found with the techniques 
described in Sections 0] and EJ The tabulated constellations have some of the best diversity 
sums and diversity products published so far. All the constellations searched by simulated 
annealing (SA) were based on the A k B k structure. For the constellations with L elements 
and parameters x, y, z being multiples of 2ir/L, they are found by geometrical methods using 
the parameterization (J4.1|) . For the constellations with L elements and parameters x,y,z 
being decimals, they are found by Brute Force with step size 0.1000 based on the same 
parameterization (j4.1lh 
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Table 4. Diversity product of 2 dimensional constellation based on weak group structure: 



Number of 
elements 


Diversity Product 


Codes and Comments 


2 


1 


x = 7T,y = 7T, z = (optimal) 


3 


V3/2 


x = 2tt/6, y = 27T/3, z = (optimal) 


4 


0.7831 


r\ c r\r\r\r\ a Ar\r\r\ 

x = 0.6000, y = 6.0000, z = 4.4000 


5 


a/5/8 


x = 2n/5, y = 87r/5, z = 4n/5 (optimal) 


8 


0.7071 


x = 2.3562, y = 3.9270, z = 4.7124 


9 


0.6524 


SA searched code 


10 


0.6124 


x = 2n/5,y = 8n/5,z = n/5 


16 


y/2/2 


x = tt/4, y = 5tt/4, z = 13tt/8 


17 


0.5255 


SA searched code 


18 


0.5207 


SA searched code 


19 


0.5128 


SA searched code 


20 


0.5011 


x = 1.6500, y = 3.7500, z = 4.0500 


24 


0.5000 


x = tt/12, y = 5tt/12, z = n/2 


37 


0.4461 


x = 2tt/37, y = 6tt/37, z = 12tt/37 


39 


0.3984 


x = 8tt/39, y = 34vr/39, z = 36tt/39 


40 


0.3931 


x = 37r/10, y = Htt/10,2 = 3tt/4 


55 


0.3874 


x = 2tt/55, y = 68tt/55, z = fa /ll 


57 


0.3764 


x = 2tt/57, y = 40tt/57, z = 48vr/57 


75 


0.3535 


x = 2tt/75, y = 98tt/75, z = 96vr/75 


85 


0.3497 


x = 267r/85, y = 94tt/85, z = 18tt/17 


91 


0.3451 


x = 2tt/91, y = 1287T/91, ^ = 42tt/91 


96 


0.3192 


x = 7tt/16, y = 29tt/16, z = tt/6 


105 


0.3116 


x = 2tt/105, y = 68tt/105, z = 84tt/105 


120 


0.3090 


x = vr/30, y = llvr/30, z = tt/4 


135 


0.2869 


x = 2tt/135, y = 28tt/135, z = 68tt/135 


145 


0.2841 


x = 2tt/145, y = 64tt/145, z = 76tt/145 


165 


0.2783 


x = 2tt/33, y = 20tt/33, z = 2tt/5 


203 


0.2603 


x = 2tt/203, y = 290tt/203, z = 70vr/203 


225 


0.2499 


x = 82vr/225,y = 118vr/225,z = 126vr/225 


217 


0.2511 


x = 2ir/2l7,y = 250tt/217,2 = 168tt/217 


225 


0.2499 


x = 827r/225, y = 1187r/225,2 = 126tt/225 


240 


0.2239 


x = 7r/40, y = 97r/40, z = 7r/6 


273 


0.2152 


x = 2tt/273, y = 208tt/273, z = 142tt/273 


295 


0.2237 


x = 14tt/295, y = 104tt/295, z = 22tt/59 


297 


0.1910 


x = 242tt/297, y = 548tt/297, z = 54tt/297 


299 


0.1858 


x = 8tt/299, y = 220tt/299, z = 18vr/299 


300 


0.1736 


x = tt/150, y = 5br/150, z = 5tt/6 
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Table 5. Diversity sum of 2 dimensional constellation based on weak group structure 



number of 
elements 


Diversity Sum 


Codes and Comments 


o 
I 


1 


x = 7r,y = tt, z = (optimal) 


3 


V3/2 


x = 271/6, y = 2,71/6, z = (optimal) 


5 


y/5/8 


x = 2tt/5, y = 8rr/5, z = 47r/5 (optimal) 


9 


3/4 


x = 107r/9, y = 47r/3, z = Ati/9 (optimal) 


16 


V2/2 


x = 7r/4, y = 57r/4, z = 137r/8 (optimal) 


18 


0.6614 


x = 4vr/9, y = 2vr/3, z = 7vr/9 


19 


0.6391 


SA searched code 


20 


0.6338 


SA searched code 


21 


0.6307 


SA searched code 


22 


0.6154 


SA searched code 


24 


0.6124 


x = n/6,y = 7r/4, z = hit/ 12 


28 


0.5996 


x = 37r/8, y = 7r/2, z = 2n/7 


30 


0.5934 


x = 47t/15, y = 7r/3, z = 7n/lh 


31 


0.5739 


SA searched code 


32 


0.5734 


SA searched code 


39 


0.5726 


x = 14tt/39, y = 40tt/39, z = 18tt/39 


40 


0.5499 


x = 3vr/20, y = 7tt/20, z = 3tt/10 


42 


0.5371 


x = An/7, y = 13tt/21, z = tt/3 


45 


0.5342 


x = 2vr/9, y = Att/9, z = Un/15 


52 


0.5332 


x = vr/13, y = 2vr/13, z = 9tt/26 


57 


0.5053 


x = 47r/57, y = 8tt/57, z = 40vr/57 


60 


0.5000 


x = vr/15, y = 4tt/15, z = 3tt/10 


64 


0.4852 


x = 3tt/16, y = 53tt/32, z = 55tt/32 


75 


0.4850 


x = 32tt/75, y = 14tt/75, z = 2vr/75 


76 


0.4672 


x = 3tt/19, y = 4vr/19, z = llvr/38 


77 


0.4595 


x = 52?r/77, y = 82tt/77, z = 60vr/77 


85 


0.4540 


x = 2tt/17, y = 8tt/17, z = 14tt/85 


87 


0.4460 


x = 52tt/87, y = 98vr/87, z = 82vr/87 


95 


0.4418 


x = 6tt/19, y = 2tt/95, z = 36vr/95 


96 


0.4390 


a; = 39tt/48, y = 5vr/12, z = llvr/24 


99 


0.4297 


x = 627r/99,y = 192tt/99,2 = 142vr/99 


105 


0.4295 


x = 27r/105,y = 16tt/105,2 = 28tt/105 


106 


0.4161 


x = 2tt/53, y = 13tt/53, z = 12tt/53 


120 


0.4156 


x = 7r/10,y = 7r/6,2; = 57r/4 


123 


0.4077 


x = 1887r/123,y = 38tt/123,2 = 182tt/123 
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number of 
elements 


Diversity Sum 


Codes and Comments 


130 


0.4071 


x = 26tt/65, y = 5tt/13, z = 2tt/13 


133 


0.3971 


x = 2tt/133, y = 2127r/133, z = 206tt/133 


138 


0.3963 


x = lQn/Q9,y = 19n/Q9,z = 4vr/69 


145 


0.3949 


x = 138tt/145, y = 22vr/145, z = 40vr/29 


148 


0.3840 


x = 5vr/74, y = 13tt/37, z = 2vr/37 


150 


0.3758 


x = tt/15, y = 8tt/75, z = 19tt/75 


155 


0.3828 


x = 2vr/5, y = 26vr/31, z = 58tt/31 


156 


0.3824 


x = 5tt/39, y = 8tt/39, z = 15tt/78 


158 


0.3823 


x = 58tt/79, y = 81tt/79, z = 64vr/79 


159 


0.3814 


x = 8tt/159, y = 64tt/159, z = 30tt/159 


160 


0.3802 


x = 69tt/80, y = 59tt/80, z = 37tt/20 


162 


0.3770 


x = 53tt/21, j/ = IOtt/9,2 = 19rr/81 


165 


0.3760 


x = 24vr/165, y = 26tt/165, z = 34tt/165 


166 


0.3699 


x = 14tt/83, y = 21tt/83, z = IOtt/83 


169 


0.3696 


x = 56tt/169, y = 76tt/169, z = 284tt/169 


171 


0.3678 


x = 32tt/171, y = 294tt/171, z = 6tt/171 


178 


0.3664 


x = 145tt/89, y = 26vr/89, z = IOtt/89 


180 


0.3636 


x = tt/9, y = 97tt/90, z = 127vr/90 


193 


0.3598 


x = 90tt/193, y = 98tt/193, z = 26tt/193 


204 


0.3566 


x = 13tt/51, y = An/51, z = 5tt/34 


208 


0.3501 


x = tt/13, y = 8tt/13, z = 65tt/104 


214 


0.3476 


x = 98tt/107, y = 67tt/107, z = 59tt/107 


220 


0.3459 


x = 197r/ll,y= 1637r/110,2 = 121tt/110 


222 


0.3438 


x = 197r/lll,y = 227r/lll,z = 157r/lll 


225 


0.3420 


x = 2tt/225, y = 52tt/225, z = 414tt/225 


234 


0.3410 


x = Air/U7,y = 24tt/117,z = 43vr/117 


240 


0.3371 


x = 7l7r/120,y = ll7r/10,2= 187tt/120 


244 


0.3335 


x = 39tt/122, y = 14tt/61, z = 20tt/61 


245 


0.3305 


x = 16tt/245, y = 186vr/245, z = 46vr/245 


248 


0.3291 


x = 103tt/124, y = 39tt/31, z = 179tt/124 


259 


0.3288 


x = 30tt/259, y = 44tt/259, z = 42vr/259 


262 


0.3274 


x = 142tt/131, y = 215tt/131, z = 87tt/131 


264 


0.3247 


x = 79tt/66, y = 129tt/66, z = 215tt/132 


276 


0.3237 


x = 23tt/138, y = 15tt/69, z = 6tt/69 


287 


0.3188 


x = 6tt/287, y = 76tt/287, z = 28tt/287 


292 


0.3164 


a; = 65tt/146, y = 14vr/73, z = 82tt/73 


295 


0.3147 


x = tt/5, y = 50tt/59, z = 22tt/59 


300 


0.3126 


x = vr/75, y = 17tt/150, z = 9tt/25 
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5.5 General Form Constellation Numerical Design 



The connection between the complex Stiefel manifold and U(M) (see the beginning of this 
section) makes clear that the techniques used above for square unitary constellations can 
be applied to design general form unitary constellations too. For simplicity we describe the 
idea with assumption T = 2M and consider the following structure: 

{A k B\A e U(T), B = ( ^ ) , k = 0, 1, ■ • • , L - 1}. 

One can check at most 2L — 1 distance calculations are needed to derive the diversity product 
(sum or function) with this algebraic structure. 

Table 6. The following tables show the constellations (M=2) found using SA. More results 
can be found in I1UI. 



size 


3 


4 


5 


6 


7 


8 


9 


diversity sum 


0.8654 


0.7901 


0.7889 


0.7652 


0.7514 


0.7422 


0.7369 


diversity product 


0.8582 


0.7424 


0.7330 


0.6450 


0.6361 


0.6216 


0.5822 



6 Fast Decoding of the Structured Constellation 

The complexity of ML decoding for unitary space time constellations increases exponentially 
with the number of antennas or the transmission rate. This will preclude its practical 
use for high transmission rates or for large number of antennas. Basically our structured 
constellations can convert the ML decoding to lattice decoding naturally, consequently they 
admit fast decoding algorithms. 

The principle of sphere decoding [S] is as follows: instead of doing an exhaustive search 
over all the lattice points, one can limit its search area to a sphere with given radius \/~C 
centered at received point. One can check the complexity of this approach in jB] and in |16j . 

We will use the A k B l structure to describe how one can apply sphere decoding algorithm 
for the demodulation based on our constellations. Suppose A has Schur decomposition A = 
Udiag (e iai ,e ia ' 2 , ■■■ , e iaM )U*, similarly assume B = Ediag (e i/3 \ e ifti , • ■ ■ , e if3M )B*. Consider 
unitary differential modulation ^H] and denote with X T the received signal at time block r. 
The ML demodulation algorithm involves the following minimization problem: 

(k, 1) = argmin \\X T - A k B l X T _ 1 \\ F . 

Algebraically one can check that 

|| X T — A k B^ 1 X T _\ \\p = \\A k X T — -B'vY" T _;jJ|p 

= ||£/diag (e- ika \ e~ ika \ • ■ • , e~ ikaM )U*X T - Vdiag (e~ m , e~ m , ■ • ■ , e- ill3M )V*X T ^\\ F 

So every entry of X T — A k B l X T ^\ is a linear combination of trigonometric functions cos or sin 
in the variables k, I, which can be viewed as lattice points. As demostrated in [20] and |16j . 
the whole demodulation task has been converted to least-squares problem. Consequently 
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our structured constellation will admit sphere decoding algorithm. In [20] a detailed study 
of the sphere decoding algorithm applied to constellations from Sp(2) was undertaken. 

The complexity (either upper bound or average complexity) of sphere decoding will de- 
pend on the dimension of the lattice. This will make the weak group structure A k B k more 
remarkable, because in this case the algorithm requires considering finding the closest point 
in a one dimensional lattice, which is very simple. 

In jS] a very interesting fast demodulation approach is proposed for diagonal space time 
constellations. The authors use numerical approximation and LLL basis reduction technique 
to reduce the decoding complexity. Note that a constellation with the weak group structure 
A k essentially is a diagonal constellation (straightforward Schur decomposition will show 
this), therefore the same technique can be applied to this structure. Most importantly 
some other algebraic structure can employ the techniques too. For instance, consider the 
A k B l C m structure. If we let I go over a large interval and let k, m stay within a small 
interval, the structure will become "almost" diagonal. For efficient decoding, one only has 
to do exhaustive search for k, m and apply the techniques for diagonal constellations to 
decode I. Although the decoding complexity will increase a little, our experiments show the 
performance will output the diagonal one remarkably. Exactly the same "almost" diagonal 
idea can be applied to other proposed structures. 

7 Conclusions and Future Work 

In this paper, we study the limiting behavior of the diversity function by either letting the 
SNR go to infinity or to zero. Respectively the diversity product and the diversity sum for 
unitary constellations are studied from the analysis of the limiting behavior. We propose 
algebraic structures, which are suitable for constructing unitary space time constellation and 
feature fast decoding algorithms. Based on the presented structure we construct unitary 
constellations using geometrical symmetry and numerical methods. For 2 dimension most 
of our codes are better or equal to the currently existing ones. For higher dimensions many 
codes with excellent diversity are found, which were never found before. Combined with 
the proposed algebraic structure the numerical methods can also be employed to optimize 
the diversity function at a certain SNR. Future work may involve analyzing the geometric 
aspects (such as geodesies, gradients and Hessians of the functions, etc) on U(M) or the 
complex Stiefel manifold. Using the optimization techniques on Riemmannian manifold to 
optimize the distance spectrum of a unitary constellation to further search good-performing 
constellations is under close investigation too. 
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